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Abstract 

Synchronization control in arrays of identical output-coupled continuous- 
time linear systems is studied. Sufficiency of new conditions for the ex- 
istence of a synchronizing feedback law are analyzed. It is shown that 
for neutrally stable systems that are detectable form their outputs, a 
linear feedback law exists under which any number of coupled systems 
synchronize provided that the (directed, weighted) graph describing the 
interconnection is fixed and connected. An algorithm generating one such 
feedback law is presented. It is also shown that for critically unstable sys- 
tems detectability is not sufficient, whereas full-state coupling is, for the 
existence of a linear feedback law that is synchronizing for all connected 
coupling configurations. 

1 Introduction 

Synchronization is, at one hand, a desired behaviour in many dynamical systems 
related to numerous technological applications (TnHHlin]; and, at the other, a 
frequently-encountered phenomenon in biology [23l[28l[l]. Aside from that, it 
is an important system theoretical topic on its own right. For instance, there is 
nothing keeping us from seeing a simple Luenberger observer [15] with decaying 
error dynamics as a two-agent system where the agents globally synchronize 
[m Ch. 5]. Due to abundance of motivation to investigate synchronization, a 
wealth of literature has been formed, mostly in recent times T, Sec. 5], [241129] . 

An essential problem in control theory is to find general conditions that im- 
ply synchronization of a number of coupled individual systems. This problem, 
which is usually studied under the name synchronization stability, has been at- 
tacked by many and from various angles. Two cases, which partly overlap, are of 
particular interest: (i) where the dynamics of individual systems are relatively 
primitive (such as that of an integrator) yet the coupling between them satis- 
fies but minimum regularity stipulations (e.g. time-varying, delayed, directed, 
need not be connected at all times); and (ii) where the individual systems are 
let be more complex, but the structure of their interaction resides in a more 
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stringent set (e.g. fixed, balanced, symmetric, with coupling strength greater 
than some threshold.) The studies concentrated on the first case have resulted 
in the emergence of the area now known as consensus in multi-agent systems 
[Tfl [TBI [HI [131 [201 [21 M US] where fairly weak conditions on the interconnection 
have been established under which the states of individual systems converge to 
a common point fixed in space. The second case has also accommodated im- 
portant theoretical developments especially by the use of tools from algebraic 
graph theory [3D]. Using Lyapunov functions, it has been shown that spectrum 
of the coupling matrix plays a crucial role in determining the stability of syn- 
chronization |331 118] notwithstanding it need not be explicitly known |5j . It 
has also been shown that passivity theory can be useful in studying stability 
provided that the coupling is symmetric [H [191 [22] • 

As alluded to earlier, Luenberger observer makes a particular example for 
synchronization. The well-known dynamics read 

X — Ax 

X ~ Ax + L{y-'Cx) 

where y = Cx is the output of the observed system x — Ax and x is the estimate 
of the actual state x. Taking observer gain as a design parameter, detectability 
of the pair (C, A) is necessary and sufficient for the existence of a linear feedback 
law L that ensures \x{t) — x{t)\ ^ as i — > oo. Once the above dynamics are 
rewritten as 
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a natural generalization of the (observer design) problem becomes apparent. 
Namely, for the following system 



Xl 



{I(S)A + T(S)LC) 



Xl 



investigate conditions on (coupling) matrix T and pair (C, A) under which one 



can find an L (synchronizing feedback law) yielding \xi{t) 



for all 



z, j (synchronization.) Most applications related to this question render the 
assumption that coupling matrix is known too restrictive. Even the size of T 
(hence the number of coupled systems) may sometimes be unknown. Therefore 
when designing a synchronizing feedback law, one is forced to not count on the 
exact knowledge of coupling matrix (interconnection.) However, something has 
to be assumed and how loose that assumption can get is yet a partly open prob- 
lem. To be precise, given a set of possible interconnections, where each element 
of the set determines not only the coupling between the systems but also the 
size of the array, a fundamental question is the following. What condition on 
pair (C, A) guarantees the existence of an L under which the systems synchro- 
nize for all coupling configurations that belong to the given set? Clearly, the 
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choice of the set shapes the possible answers to this question. When the set is 
the collection of all interconnections with coupling strength greater than some 
(a priori known) positive threshold, the answer (which can be extracted, for in- 
stance, from [32]) is as simple as one can expect: that (C, A) is detectable, just 
as is the case with Luenberger observer. Unlike Luenberger observer though, 
mere that A — LC is Hurwitz is not sufhcient for synchronization. In attempt to 
make the general picture closer to complete (see Theorem [Ij, in this paper we 
focus on the set of all interconnections with connected graphs (no assumption 
on the strength of coupling nor on the symmetry or balancedness of the graph) 
and study new conditions on (C, A) which do (not) guarantee the existence of 
a synchronizing feedback law with respect to that set. The main contribution 
of this paper is in establishing the following results: 

(b) If (C, A) is detectable with A neutrally stable then there exists a synchro- 
nizing linear feedback law. (We also provide an algorithm to explicitly 
compute one such feedback law.) 

(c) If A is critically unstable (e.g. double integrator) and C is full column 
rank then there exists a synchronizing linear feedback law. 

(f) If (C, A) is detectable with A critically unstable then, in general, there 
does not exist a synchronizing linear feedback law. 

For the discrete-time version of (b) see [26] . 

In the remainder of the paper we first provide basic notation and definitions. 
In Section[3]we define synchronizability with respect to a set of interconnections 
and formalize the problem through that definition. Following the problem state- 
ment we present the main theorem as a 8-item list of sufficient and nonsufficient 
conditions for synchronizability. After stating the main theorem (Theorem [1} 
we set out to prove what it claims. The sufficiency statements of Theorem [T] are 
demonstrated in Section |4] where we also show that coupled harmonic oscillators 
synchronize if the coupling is via a connected graph (Corollary[T]) To the best of 
our knowledge, establishing synchronization of harmonic oscillators without any 
symmetry, balancedness, or strong coupling assumption on the interconnection 
is new. Section [5] includes the proofs of nonsufficiency statements of Theoremll] 
which are given in order to give a measure on how tight (close to necessary) the 
sufiicient conditions are. We spend a few words on a dual problem in Section [6] 
after which we conclude. 



2 Notation and definitions 

Let N denote the set of nonnegative integers and M>o set of nonnegative real 
numbers. Let | • | denote 2-norm. For A £ C let Re(A) denote the real part of A. 
Identity matrix in R"^" is denoted by /„ and zero matrix in R™^" by 0„ixn- 
Conjugate transpose of a matrix A is denoted by A^ . A matrix A S C"^" 
is Hurwitz if all of its eigenvalues have strictly negative real parts. A matrix 
S e M"""" is skew- symmetric if 5* + S*^ = 0. For A e R"""" and C G R™^", pair 
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(C, A) is detectable (in the continuous-time sense) if that Ce 



for some 



x G R" and for alH > imphes hmj^oo e'^^x = 0. Matrix A e R"^" is neutrally 
stable (in the continuous-time sense) if it has no eigenvalue with positive real 
part and the Jordan block corresponding to any eigenvalue on the imaginary 
axis is of size one. Let 1 € R'' denote the vector with all entries equal to one. 
Kronecker product of A e R™><" and B eW^'^ is 



A(RiB 



aiiB 



O-mlB 



ainB 



rB 



Kronecker product comes with the properties {A ^ B){C ^ D) ~ (AC) (BD) 
(provided that products AC and BD are allowed) Aig) B + Aig) C ^ A® {B + C) 
(for B and C that are of same size) and {A (g) B)^ — A^ (g) B^ . 

A (directed) graph is a pair [Af, £) where A/" is a nonempty finite set (of nodes) 
and £ is a finite collection of ordered pairs [edges) (n^, Uj) with ni, Uj € £. A 
path from ni to ni is a sequence of nodes (ni, n2, . . . , ng) such that (n^, n^+i) 
is an edge for i e {1, 2, 1}. A graph is connected if it has a node to 

which there exists a path from every other nodel3 Fig. [1] illustrates two graphs, 
where one is connected and the other is not. 




n2 




ni 



Go 



Figure 1: Schematic representations of two graphs Gi = (Af, £i) and G2 
(A/", £2) where A/" = {ni, 712, na}, = {(ni, 712), (^2, rts), ("3, "2)} and £2 
{(712, ni), (n2, 713)}. Gi is connected, G2 is not. 



A matrix F := [7^] g R^^^ describes (is) an interconnection if 7^^ > for 
i ^ j and = — ^j^i lij ■ It immediately follows that A = is an eigenvalue 
with eigenvector 1 (i.e. Fl = 0.) The graph of F is the pair [Af, £) where 
Af = {ni, 772, • ■ • J np} and (71,;, nj) S f iff 7^ > 0. Interconnection F is said to 
be connected if its graph is connected. 

For connected F, eigenvalue A = is distinct and all the other eigenvalues 
have real parts strictly negative. When we write Re(A2(F)) we mean the real 
part of a nonzero eigenvalue of F closest to the imaginary axis. Let r'^ be the 

^Note that this definition of connectedness for directed graphs is weaker than strong con- 
nectivity and stronger than weak connectivity. 
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left eigenvector of eigenvalue A = (i.e. r^T = 0) with r^l = 1. Then is 
unique and satisfies limt^oo e'"* = Ir'^ . 

Given maps : R>o — > M" for i = 1, 2, . . . , p and a map ^ : M>o — » K", 
the elements of the set {^i(-) : i — 1, 2, . . . , p} are said to synchronize to ^(•) 
if \^i{t) — ^{t)\ ^ as t — + 00 for all i. The elements of the set {^i(-) : i = 
1, 2, . . . , are said to synchronize if they synchronize to some ^(•). 

Let V denote the set of all pairs (C, A) where matrix C and square matrix 
A, both real, have the same number of columns. We define the following subsets 



• .Ah: set of all pairs (C, A) with A Hurwitz. 

• A-^: set of all pairs (C, A) with A neutrally stable. 

• Aj: set of all pairs (C, A) with A having no eigenvalue with positive real 



• Of' set of all pairs (C, A) with C full column rank. 

• Op: set of all detectable pairs. 

Few remarks are in order regarding the above definitions. Note that C 
Af^ C Aj and Op C Op. Set Aj allows A matrices with Jordan blocks of 
arbitrary size with eigenvalues on the imaginary axis. (Hence the subscript J.) 
Therefore there are pairs (C, A) in Aj with critically unstable A matrices. For 
a pair (C, ^4) e Op if C is an output matrix of some system of order n then 
there exists matrix L such that LC = /„, i.e. full state information is instantly 
available at the output. (Hence the subscript F.) For an arbitrary (C, A) G Op, 
however, the state information is only partially available at the output. (Hence 
the subscript P.) We also need the following definitions. 

• Q>s- set of all connected interconnections with |Re(A2(r))| > 5 > 0. 

• Q>o- set of all connected interconnections. 

• Q>o- set of all interconnections. 
Note that g>s C ^>o C g>o. 

3 Problem statement and main theorem 

For a given interconnection F e ffi^^^, let an array of p linear systems be 



ofV. 



part. 



Vi = Cxi 




(la) 
(lb) 
(Ic) 
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where Xi G is the state, Ui G R" is the input, yi G M™ is the output, and 
Zi G K™ is the coupling of the ith system ioi i — 1, 2, . . . , p. Matrices A and C 
are of proper dimensions. The solution of ith system at time t > is denoted 
by Xi{t). 

Definition 1 Given A gM.^^"^, C G M™'*", and set of interconnections S; pair 
(C, A) is said to be synchronizable with respect to S if there exists a linear 
feedback law L G M"^™ such that for each T G S solutions Xi{-) of array (fT|) 
with Ui — Lzi synchronize for all initial conditions. 

When is pair (C, A) synchronizable with respect to interconnection set QyQ ? 

The above question is what we are mainly concerned with in this paper. 
Our objective is to find sufficient conditions on (C, A) for synchronizability with 
respect to set of all connected interconnections and examine their slackness. To 
this end, we adopt a constructive approach, i.e. we compute a synchronizing 
feedback law L whenever possible. We also provide the explicit trajectory that 
the solutions synchronize to (if they do synchronize.) By adding our findings 
to already known results on synchronization of coupled linear systems we aim 
to reach a clearer picture of the problem with fewer missing pieces, which is 
recapitulated in the below theorem. 

Theorem 1 Let us be given 5 > Q. We have the following. 

(a) All (C, A) G An are synchronizable with respect to t/>o- 

(b) All (C, A) G .4n n Op are synchronizable with respect to Q^q. 

(c) All (C, A) G A3 n Op o.i'e synchronizable with respect to QyQ. 

(d) All (C, A) G Op are synchronizable with respect to Q>s- 

(e) Not all (C, A) G .4n H Op are synchronizable with respect to Q>o. 

(f) Not all (C, A) G Aj n Op are synchronizable with respect to t?>o- 

(g) Not all (C, A) G Op are synchronizable with respect to Q^q. 

(h) Not all (C, A) G V are synchronizable with respect to G>s. 

As mentioned earlier, statement (d) has been known for a while. General 
results pertaining to that case go as far back as [33] . Among the remaining 
statements, (a), (e), (g), and (h) are straightforward if not obvious. Still, for 
the sake of completeness, we establish each of the eight statements of Theorem[T] 
by a series of lemmas. 
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4 Cases when synchronizing feedback exists 



This section is dedicated to prove the existence statements of Theorem [1] which 
are analysis results. However, what is in here is more than mere showing the 
existence of some synchronizing feedback law L under different assumptions. 
For each case we actually propose/compute an L. We also explicitly provide 
the trajectory that systems will synchronize to under the proposed feedback 
law. Therefore what will have been solved at the end is a synthesis problem. 

4.1 System matrix Hurwitz 

Let us be given some interconnection T. Consider ([1]) under the feedback con- 
nection Ui — Lzi for some L. When system matrices A are Hurwitz, the trivial 
choice L = decouples the systems and (t) — > as i ^ oo for all initial condi- 
tions. Therefore synchronization accrues. We formalize this trivial observation 
to the following result and, consequently, establish Theorem [TJa). 

Lemma 1 Given (C, A) e Ab_, let L := 0. Then for all T G Q>q solutions Xi{-) 
of systems ([T|) with Ui — Lzi synchronize to x{t) = 0. 

4.2 System matrix neutrally stable 

The second statement of Theorem [T] is concerned with neutrally stable system 
matrices. Here we choose to first study a subcase (Proposition [T]) where the 
system matrices are skew-symmetric and then generalize what we obtain. For 
some interconnection F — [7^] consider the following coupled systems 

p 

i^ = Si^+H^HY,^,,{^,-i{), i = l, 2, (2) 

where £,i e M" is the state of the iih system, S £ M""=", and H e M™^". We 
make the following assumptions on systems ([2]) which will henceforth hold. 
(Al) S is skew-symmetric. 
(A2) (H, S) is observable. 
(A3) F is connected. 

Proposition 1 Consider systems Solutions ^i(-) synchronize to 

ep(o). 

where r G W is such that r'^T — and r"^! — 1. 
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Proof. Consider matrix F — Ir'^. Observe that (T — Ir'^)'' = T'^ + (-l)'^lr'^ 
for fc € N. For i € R therefore we can write 

e^^-i'-")* = /, + f(F-lr^) + ^(F-lr^)2 + ... 

^ (/, + .F4F^ + ...)-(n.-4l.- + ...) 

= e^' ~ {1 - e-')lr^ . 

Consequently hmt^oo e^'"^-'-'' ^* = 0. We deduce therefore that F — Ir^ is Hur- 
witz. Since F— Ir'^ is Hurwitz, there exist symmetric positive definite matrices 
P,Qe RP^'P such that 

-Q={r^lr^fP + P{T-lr'^). (3) 

Define positive semidefinite matrices P := {Ip — Ir^)^ P{Ip — Ir^) and Q := 
{Ip — lr'^y Q{Ip — lr'^). Now pre- and post-multiply equation ([3]) by (/p — Ir^)^ 
and {Ip — Ir^), respectively. We obtain 

-Q = {Ip-lr^f{T-lr^fP{Ip-lr^) 

+ [I^-lr^fp{T-lr^){Ip-lr^) 
= F^P(/p - Ir^) + {Ip - lr^)^PF 

= F^(/p - lr'^ fP{Ip - Ir^) + {Ip - lr^ fP{Ip - lr^)F 
= F^P + PF . 

We now stack the individual system states to obtain x := [^^ C J ■ • • ^pV ■ We 
can then cast ([2]) into 

i = (/p ® 5 + F ® H^H)Ti . (4) 

Define V -.W"^ ^ ]R>o as ^(x) x^(P«)/„)x. Differentiating y(x(i)) with 
respect to time we obtain 

y(x) = x^{Ip(g>S^ + (g>H'^H){P(g)In)^ 

-hx^(P (8) /„)(/p ® + F «) i/^i/)x 
= x^(P(»(S'^ + 5) + (F^P + PF)®i7^i7)x 
= -x^{Q(g> H^H)^. (5) 

Thence ^(x) < for both Q and H^H (and consequently their Kronecker 
product) are positive semidefinite. 

Given some ( G MP^, let X C W"^ be the closure of the set of all points rj 
such that T] = {Ir^ e'^*)C for some t > 0. Set X is compact for it is closed by 
definition and bounded due to that C is fixed and S" is a neutrally-stable matrix. 
Having defined X, we now define 

n:={Tje RP" : (Ir^ ® In)v & X , V{r,) < V{0} ■ 
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Let us show that f2 is forward invariant. Observe that 

^ ((Ir^ (g) InMt)) = (Ir^ «) /„) (/p (g) + r (g) H^H)x{t) 

= (Ir^ g) S" + Ir'^r (g) i/^i/)x(t) 

= (Ir^ (g 5)x(i) 

= (/p(gS')(lr'^(g/„)x(t). 

We therefore have 

{Ir^ g) /„)x(i) = (Ir^ g) e^*)x(0) (6) 

which in turn imphes that if (Ir^ g) /„)x(0) G A" then (Ir-^ g) /„)x(t) e A" for 
all f > 0. Likewise, if y(x(0)) < V{C) then y(x(t)) < 1/(C) for aU t > thanks 
to ([5]). As a result, if x(0) € fl then x(i) G for all t > 0, that is, SI is forward 
invariant with respect to (g]). 

Set n is closed by construction. To show that it is compact therefore all we 
need to do is to establish its boundedness. Let 

a := sup |77 - (Ir^ (g /„)ry| . 

v{v)<viC) 

If we go back to the definition of V wc immediately see that a < oo. Now let 

b :— sup \uj\ . 

Since X is bounded, 6 < oo as well. Now, given any 77 e O we have |?7 — (Ir'^ (g 
In)'>l\ ^ Hence we can write 

\v\ < a+|(lr^g)/„)?7| 
< a + sup |a;| 

= a + b. 

Therefore f2 is bounded. Having shown that is forward invariant and compact, 
we can now invoke LaSalle's invariance principle il4, Thm. 3.4] and claim that 
any solution starting in approaches to the largest invariant set W C G SI : 

viv) = 0}. 

Let now 7y(-) be a solution of Q such that ri{t) G W for all t > 0. Given 
some T > 0, since Vijiir)) = 0, we can write 

which implies, since Q is positive definite, that either ((/p — Ir-^) ® In)rj{T) ~ 
or {Ip g) H)rj{T) = 0. Suppose now that 

((/p-lr^)®/„MT) ^0. (7) 
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Continuity of implies that there exists 6 > such that {{Ip — lr'^)CSiIn)r]{t) ^ 
for t G [t, t + S]. Therefore we must have (/p (g) H)T]{t) = for i e [r, r + J]. 
However, observabihty of pair {H, S) stipulates that r]{t) — for t £ [t, t + 5] 
which contradicts We then deduce ((/p - Ir^) (g) = for all t > 0. 

Therefore W C {lo e n : u = (Ir'^ J„)w} = X. 

Let us now be given any solution x(-) of dl]). Since C that we used to 
construct fl was arbitrary, without loss of generality, we can take x(0) = C- 
That x(0) € implies that x(t) approaches X a.s t ^ oo. Therefore we are 
allowed to write 

= lim (x(t) - {Ir'^ ® InMt)) 
= lim (x(t) - (Ir"^ (E) e^*)x(0)) 

where we used (O- ■ 
The following result, which establishes synchronization of coupled harmonic 
oscillators, comes as a byproduct of Proposition [T] 

Corollary 1 Consider p coupled harmonic oscillators (inM.'^) described by 

— yi 

p 

Oscillators synchronize for all connected interconnections [7^] . 

We will use the next fact in Algorithm [1] to construct a feedback law for syn- 
chronization. 

Fact 1 Let F G R"^" be a neutrally-stable matrix with all its eigenvalues re- 
siding on the imaginary axis. Then 

P := lim / e^^^e^^dT (8) 

is well-defined and symmetric positive definite. It also satisfies PF + F^P = 0. 

Proof. Matrix F is similar to a skew-symmetric matrix. Therefore e^* is 
(almost) periodic [57j. Periodicity directly yields that limit in ^ exists, that 
is, P is well-defined. Similarity to a skew-symmetric matrix also brings that 
inftgR |e^*| > and sup^gjj |e^*| < 00. Same goes for F^. Therefore there exist 
scalars a, b > such that a/„ < *e^* < &/„ for all t G M. We can then write 

ain < / e^^^e^^dT < bin 
Jo 



10 



for all t > 0. Therefore P is positive definite. Symmetricity of P comes by 
construction. Finally, observe that 



\PF + F' P\ 



= lim t^^ 

>oo 

= lim 

t—^OO 



F + F^e 



dT 



-( 



< lim t 

t—*oo 



< lim t-^{b + l) 



whence the result follows. 



Algorithm 1 Given A e K"^" that is neutrally stable and C e R™^", we 
obtain L S ^5 follows. Let ui < n be the number of eigenvalues of A that 

reside on the imaginary axis. Let 712 := n — ui. Lf ui = Q, then let L :— 0; else 
construct L according to the following steps. 
Step 1: Choose U € and W G R"''"^ satisfying 



[U W]-^A[U W] 



F 
G 



where all the eigenvalues of F ^ ]g^"ix»ii kave zero real parts. 
Step 2: Obtain P G R"i><"i from F by 
Step 3: Finally let L := UP-^{CU)'^. 

Below result establishes Theorem [TJb). 

Lemma 2 Given (C, A) G H Op, let L be constructed according to Algo- 
rithm]^ Then for all T G Q>o solutions Xi{-) of systems ([T]) with m — Lzi 
synchronize to 



x{t) := (r^ 



where r G R^ is such that r^^T — and r"^! — 1. 

Proof. Let the variables that are not introduced here be defined as in Al- 
gorithm [TJ Let H CUP-^/^ and S := pi/2^p-i/2. Then {H, S) is ob- 
servable for (C, A) is detectable. Also, note that S is skew-symmetric due to 
PF + F'Fp = 0. 

We let W G R"!""" and G R"^""" be such that 

"C/f 



= [U W] 
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Note then that WU = = U^'W = 0, and W'fU = 0. Since 

Ui = Lzi, from ([1]) we obtain 



Axi + LC ^ (xj - Xi) 

3 = i 



Let now e M"! and ry, e R"^ be 



"6' 




■pl/2 " 




"c/t" 






/„, 







(9) 



(10) 



Combining ([9]) and (fTO]) we can write 

p p 

77, = (12) 
Let r e KP^P be connected and r e Rp be such that r^F = 0. Then define Ui : 



as LUi (t) 



,TlT 



^i(i) for z = 1, 2, . . . , p. Let w := [wf ^2 • ■ • ] 



and V := [yyf rj^ ■ ■ ■ rj^Y" ■ Starting from pT|) and (fT^ we can write 

w(i) = (r(g)e"^*7J^iJe^*)w(i) + (r(g)e"^*ij'^CVl/e'^*)v(0). 

Thence 



w(t) = $(t, O)w(O) + 
where 



/ $(t, t) (r ® e-^^H^CWe^'')dT 
Jo 



v(0) (13) 



$(t, t) := exp (r (g) e-^°'H^He^°') da^ 

is the state transition matrix [3] . From Proposition [T] we can deduce that 
^>(t, t) is uniformly bounded for aU t and r. Also, for any fixed r we have 
limf^oo $(i, t) = l?"^ (S" Moreover, e'^* is uniformly bounded for all t, and 
e*^* decays exponentially as t — + cxd for G is Hurwitz. Therefore we can write 

lim / <^{t,T){T ®e-^''H'^CWe^'')dT 

= [ f lim $(t, r)) (r (g> e-'^^H'^CWe^^)dT 

/>OC 

Jo 



= 0. 



Then, by p3|) . we can write 



lim w(t) = (Ir-^ (g) /„Jw(0) . 

t — >-oo 
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Therefore solutions ^i(-) synchronize to {r'^(^e^*)w{0). Moreover, hmt^oo v(t) = 
for G is Hurwitz. Hence we can say that solutions rii{-) synchronize to 
(r^ (8) e'^*)v(0). As a result, solutions Xi{-) synchronize to 

xi(0)" 
xp{0)_ 



Hence the result. ■ 



UP 



-1/2 



W 







pl/2[/t 



xi(0) 
a;p(0) 



4.3 System matrix critically unstable 

The third statement of Theorem[T]is not surprising. The naive reasoning, which 
we will later have no difhculty in formalizing, is that when the output matrix C 
is full column rank, there exists an L such that LC equals the identity matrix. 
Then, see there will be two forces driving the array. One is due to system 
matrices A, the other due to coupling between the systems. When A matrices 
have no eigenvalues with positive real part, even if the systems tend to move 
away from each other, that tendency will not be of exponential nature. How- 
ever, if r is connected, that apart-drifting behaviour will be dominated by the 
exponential attraction introduced by the coupling between the systems. The 
result will be synchronization of the solutions. Having said that, now we can 
state our main reason for including this existence statement in Theorem [TJ to 
emphasize the sixth statement of Theorem [1] which is a nonexistence result. 
When availability of full-state information is replaced with detectability, which 
in all the other cases we study causes no problem, synchronizability is lost for 
critically unstable systems. That is. Theorem [TJf) holds despite Theorem (Tfc) , 
which we find counterintuitive. Below we establish Theorem [Ijc). 

Lemma 3 Given (C, A) e Aj n Op, let L ~ {C^C)-^C^. Then for all T G 
5>o solutions Xi{-) of systems ^ with Ui — Lzi synchronize to 

xiioy 

where r E is such that r'^T — and r^l = 1. 
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Proof. For Ui — Lzi and some connected F G Rp'^p, with r G Rp such that 
r-'^F = and r^l = 1, from ([T]) one can write 

p 

ii — + LC ^ ^ 7ij Xi) 

j=i 

p 

= Axi + ^ 7y {xj - 
j=i 

since LC = In- Letting x :— [xJ' x'^ . . . x'^Y' we obtain 

x= (/p(8) A + r«)/„)x 

which impHes 

x(i) - (er*®e'4*)x(0). 

Let us define 

:= x(t)-(lr^®e'4*)x(0) 

= (e"®e^*)x(0)- (lr^®e'4*)x(0) 

= ((e^*-lr^)®e'^*)x(0). (14) 

From the proof of Theorem[2] we know that e'"* — Ir-^ = e'^'""-'-'' ^* + e~*lr^ and 
that F — Ir-^ is Hurwitz. Also recall that A does not have any eigenvalues with 
positive real part. Therefore there exist Mg, Ma, A > 0, and integer fc < n — 1 
such that je^'"^-'-'" •'*| < AlgC^-^^ and jc^*! < Mat''- Then we can proceed from 
(0 as 

i{t) = ((e^* - Ir^) (g, e^*)x(0) 

= (ef^-^"-")* ® e^*)x(0) + e-*(lr^ ® e^*)x(0) 

whence 

\m\ < • |e^*| • |x(0)| + e-*|lr^| • |e^*| • |x(0)| 

< Mge-^'Mj''\xiO)\ + e-*A4t^-|lr^| • |x(0)| 

whence limt^oo = 0. Hence the result. ■ 

4.4 System matrix arbitrary 

The fourth statement of Theorem [1] is a special case of a well-studied problem 
|331 1181 [31 1321 131j which is based on the idea that identical systems synchronize 
under strong enough coupling. For coherence we regenerate a proof here. We 
first borrow a well-known result from optimal control theory [21 : Given a de- 
tectable pair (C, A), where C G R™^" and A e R"^", the following algebraic 
Riccati equation 

AP + PA^ + /„ - PC'^CP = (15) 
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has a (unique) solution P = > 0. One can rewrite as 

{A - PC'^C)P + P{A - PC^Cf + {In + PC^CP) = 
whence we infer that A ~ PC^C is Hurwitz. 



Claim 1 LetC G and A S R"^" satisfy (flSl) for some symmetric positive 

definite P. Then for all a > 1 and uj £ R matrix A — (a+juj)PC'^C is Hurwitz. 

Proof. Let e cr - 1 > 0. Write 

{A-{a+ juj)PC'^C)P + P{A -{<T + iuj)PC^C)" 

= (A - (ct + juj)PC^C)P + P{A -{a- juj)PC'^Cf 

= (A - (1 + e)PC^C)P + P{A - (1 + e)PC^Cf 

= (A - PC^C)P + P{A - PC^Cf - 2ePC^CP 

= -In-il + 2e)PC^CP. (16) 

Observe that (jl6p is nothing but (complex) Lyapunov equation. ■ 
Below result establishes Theorem [lid). 

Lemma 4 Given S > and (C, A) e Op, let L max{l, S~^}PC'^ where P 
is the solution to (jlSp . T/ien for all T G tj>5 solutions Xi{-) of systems ([T]) with 
Ui = Lzi synchronize to 



x{t) {r' (S) e^*) 



a;i(0) 



XpiO) 



where r G is smc/i </ia< r'^T — and r-'^l = 1. 
Proof. Given T G G>5, from ([T]) we can write 

p 

Xi = Axi + ^ 7y (a^j - Xi) . 



(17) 



Stack individual system states as x 
write 



r T T 
[Xi X2 



Xp]'^ . Then from (fTTl) we can 



X = (/p (g) A + r (8) LC)x . 



(18) 



Now let Y G CP^(J'-i\ W G C(P-i)^P, G CP'^p, and upper triangular A G 



1/ = [1 r] , 
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and 



••• 
A 



Note that the diagonal entries of A are nothing but the nonzero eigenvalues of 
r which we know have real parts no greater than —6. Engage the change of 
variables v :— {V~^ (E) /n)x and modify (US]) first into 



and then into 



v^{Ip(g)A + y-^TV (g) LC)w 

Onx (p— l)n 



A 

0(p— l)nxn 



p-1 



1 A + A LC 



(19) 



Observe that Ip-i ® A + A® LC is upper block triangular with (block) diagonal 
entries of the form A + \i max{l, S^^^PC^ C for i = 2, 3, . . . , p with Re(Ai) < 
—5. Claim [T] implies therefore that Ip-i ® A + A ® LC is Hurwitz. Thus (fTI 
implies 



lim 



(p— l)n X n 



Onx (p— l)n 
0(p— l)nx (p— l)n 



v(0) 







which yields 



lim |x(t) - (Ir^ 



)x(0)| =0. 



Hence the result. 



5 Cases when no synchronizing feedback exists 

In this section the nonsufficiency statements of the main theorem are demon- 
strated. Below we establish Theorem [Ije). 

Lemma 5 There exists (C, A) G ^nHOf that is not synchronizable with respect 
to Q>Q. 

Proof. Suppose not. Then for C = 1 and A^ (note that (C, A) e A^HOf) 
there exists L G K such that for all F G G>o solutions of systems ([ij with 
Ui = Lzi synchronize for all initial conditions. However for F := 02x2 G G>o 
we have = 0, for i = 1, 2, regardless of L. Therefore solutions Xi{-) do not 
synchronize for Xi(0) ^ 2^2(0). The result follows by contradiction. ■ 
Next result establishes Theorem [IJf ) . It emphasizes the nonsufficiency of 
partial-state coupling for synchronizability of critically unstable systems with 
respect to set all connected interconnections. As shown by Lemma [3] this lost 
synchronizability can be restored under full-state coupling. 
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Lemma 6 There exists (C, A) G AjHOp that is not synchronizable with respect 
to Q^o- 



1 




Proof. Suppose not. Let 

C := [1 0] , A 
Note that (C, A) e Aj nOp. Then there exists L e M^""^ 

_P 

such that for all T e tJ>o solutions of systems ([ij with Ui = Lzi synchronize for 
all initial conditions. That is equivalent, by |18J, to that A + \LC is Hurwitz 
for all A e {7 e C : 7^ 7 is an eigenvalue of some F € G>q} C Note that 



whence 



r e ^>o =^ e g^o Vr > o 



Ae£=^rAe£ Vr>0. 



(20) 



Observe that any F € tJ>o n R^^^ has a negative (real) eigenvalue. By ([20|) 
therefore £ includes negative real line. Then a simple calculation shows that 
both £ and p need be strictly positive which lets us, thanks to (jSO]) . take = 2 
without loss of generality. Thence A + \LC reduces to the following matrix 



2A 1 
\p 



(21) 



which must be Hurwitz for all A G £. The roots of the characteristic equation 
of the matrix in (|2T|) is 



si,2 A± VA2+Ap- 
Then for 6 e (7r/2, 37r/2) we can write 



1 



R-e(si_2)|_^_gje = COS 9 ± Y COS 29 + p cos 9 + \/ (cos 29 + p cos 0)^ + (sin 26 + p sin ( 



cos I 



V2 

'± y'cos 26* + p cos 6* + Vl + + 2/3 cos 6* . 



Observe that we can choose 9 e (7r/2, tt] and m > such that the following 
inequalities are satisfied 



1 + -^T < v/l + P^ + 2/9COs6l 

—pcoa9 < — ^ 

- cos 6* < — 
m 
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for all 6 E (7r/2, 6]. Then for all 6 E {tt/2, 9] one of the roots, say si, satisfies 
Re(si) " 



cos 61 + ^\/cos26' + pcos6'+ ^/l + f^+2pcos9 



> cos ( 



^\/cos26'- + 1 + 
y/2 V m'^ TO^ 



V2 

> cos e + — 

m 

> P- 
m 

> 



C0S61 + ^i/2cos2 6i + ^ 



which implies, since no eigenvalue of matrix (|2ip can have positive real part, 
e^'^ ^ £ for 6* e (7r/2, Due to ^ therefore the wedge W := {re^'^ : 6* e 
{tt/2, 9], r > 0} and £ are disjoint, i.e. 



(22) 



For p e N>i let Tp e M^^p be 



-1 1 
-1 1 





1 



Observe for all p that we have Tp E G>o and that Xp := e-' p — 1 is an eigenvalue 
of Tp with corresponding eigenvector 

"(Ap + l)i^ 
(Ap + 1)2 

{Xp + iy 

Note that for p large enough Xp E W which contradicts (|22|) for Fj, G ^/>o- 
Hence the result. ■ 
Below we establish Theorem [IJg). 

Lemma 7 There exists (C, A) E Of that is not synchronizable with respect to 

G>Q- 

Proof. Suppose not. Then for C = 1 and A — 1 (note that (C, A) E Of) there 
exists L E M. such that for all F E G>o solutions of systems (P) with m = Lzi 
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synchronize for all initial conditions. Choose e > such that 1 — eL > 0. Then 
let 



r 







Observe that F e G>q- Consider systems ([TJ under interconnection F with 
Ui = Lzi, for i = 1, 2. We can write xi — X2 — {\ — eL){xi — X2). Therefore 
solutions Xi{-) do not synchronize unless a;i(0) = a;2(0)- Hence the result by 



Below lemma yields Theorem (rfh). 
Lemma 8 For no 6 > pair (0, 0) is synchronizable with respect to G>5- 

6 Dual problem 

For a given interconnection F e Rp^*', consider the following array of p identical 
linear systems 



where Xi G R", Ui S R"", and Zi G R" for i 1, 2, . . . , p. Matrices A and B are 
of proper dimensions. The duality between observability and controllability for 
linear systems readily yields the following theorem by which the results in this 
paper can be extended for the synchronization of the arrays of coupled linear 
systems depicted by ([23]). 

Theorem 2 Let A £ R"^", B e R"^"*, and set of interconnections S be such 
that {B^ , A'^) is synchronizable with respect to S. Then and only then there 
exists a linear feedback law K e R™^" such that for each F e 5 solutions Xi{-) 
of array ([23|) with m = Kzi synchronize for all initial conditions. 

7 Conclusion 

For arrays of identical output-coupled linear systems we have investigated the 
sufficiency of certain conditions on system matrix A and output matrix C for 
existence of a feedback law under which the systems synchronize for all coupling 
configurations with connected graphs. Namely, we have filled in the previously 
missing pieces (the boxes indicated with a question mark) of the chart given 
in Fig. [21 In addition, for each case corresponding to a box with "/" we have 
designed a synchronizing feedback law. 



contradiction. 




Axi + Bui 



(23a) 
(23b) 
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Figure 2: Sufficiency of certain conditions on (C, A) for synchronizability with 
respect to different sets of interconnections F. 
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